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Abstract — Parameter estimation in cognitive communications 
can be formulated as a multi-user estimation problem, which is 
solvable under maximum likelihood solution but involves high 
computational complexity. This paper presents a time-sharing 
and interference mitigation based EKF (Extended Kalman Filter) 
design for joint CFO (carrier frequency offset) and channel 
estimation at multiple cognitive users. The key objective is to 
realize low implementation complexity by decomposing high- 
dimensional parameters into multiple separate low-dimensional 
estimation problems, which can be solved in a time-shared 
manner via pipelining operation. We first present a basic EKF 
design that estimates the parameters from one TX user to one 
RX antenna. Then such basic design is time-shared and reused 
to estimate parameters from multiple TX users to multiple RX 
antennas. Meanwhile, we use interference mitigation module 
to cancel the co-channel interference at each RX sample. In 
addition, we further propose adaptive noise variance tracking 
module to improve the estimation performance. The proposed 
design enjoys low delay and low buffer size (because of its online 
real-time processing), as well as low implementation complexity 
(because of time-sharing and pipeling design). Its estimation 
performance is verified to be close to Cramer-Rao bound. 

I. Introduction 

Cognitive communication system is widely accepted as 
a perspective way in increasing the spectrum efficiency of 
wireless networks, where primary links and secondary links 
can usually co-exist in the network, resulting in an interference 
limited environment. Parameter estimation in cognitive com- 
munications is a challenging problem because of (i) the exis- 
tence of co-channel interference, and (ii) the high-dimensional 
parameters from multiple TX users to multiple RX anten- 
nas. In particular, note that different TX users often have 
independent carrier frequency offset (CFO) values (including 
both oscillator offsets and doppler offsets), which usually 
introduce serious nonlinear components within the observed 
signal, complicating the estimation problem. Meanwhile, chan- 
nel responses from multiple TX users to multiple RX antennas 
can result in a set of high-dimensional parameters, which are 
also difficult to estimate. Finally, due to the existence of multi- 
user interference, CFO and channel parameters usually have 
to be treated together and be estimated in a joint way so as 
to approach the optimal performance, which further increases 
the estimation complexity. 

Without loss of generality, this paper assumes Orthogonal 
Frequency Division Multiplexing (OFDM) system, which is 



an overwhelming choice for modern wireless systems. The 
classical CFO and channel estimation method in a single-user 
OFDM system is based on two repeated training symbols 
|fl~). h has low implementation complexity and near-optimal 
performance, but only applies to a single-user scenaricQ. In 
multi-user OFDMA systems with unique subcarrier set per 
TX user, CFO and channel parameters can be recovered by 
exploiting distinct subcarrier structures among TX users (|2), 
101). But this method requires separate subcarrier allocation for 
different users. Consequently, in a general multi-user cognitive 
system without specific subcarrier allocation per user and with 
overlapped training symbols, ML and EM related methods 
seem to be the only applicable choice, where all TX users' 
parameters have to be formulated into a maximum likelihood 
(ML) estimation problem 2), which is solvable under Ex- 
pectation Maximization (EM) method [5] in an iterative way. 
However, since the entire OFDM block is stored offline and 
is iteratively processed multiple times, these ML and EM 
approaches often require high computational complexity and 
high processing delay. 

Based on above considerations, this paper will focus on 
using Kalman filter structure to estimate the CFO and channel 
parameters in multi-user cognitive communications. Our major 
objective is to achieve low-complexity and low-delay estima- 
tion performance in cognitive systems. In general, Kalman 
filter is a good candidate for low delay and low complexity pa- 
rameter estimation primarily due to its real-timing processing 
property. It has been conventionally used for CFO and channel 
estimation in multi-user OFDM systems, e.g., the FFT-Block 
EKF design in |6|, the parallel EKF design in [7], and the 
particle filter design in [7|. However, these existing designs 
inherently suffer from multiple issues related to complexity, 
delay and buffer size as follows: 

1) Block EKF design in [6 1 operates on an FFT-block basis, 
which grows increasingly complex as FFT size becomes 
large (e.g., 2048 FFT size). Also, parameters estimated 
in this method are handled in a high dimension manner. 

2) Parallel EKF design in [7 1 also operates on an FFT-block 
basis, marking it complex under large FFT size. Parame- 
ters in this design are jointly estimated by calculating the 

'it is also applicable for multi-user scenario with non-overlapping training 
symbols, but this is not the case considered in this paper. 



covariance information between different users, leading 
to a high matrix dimension. 
3) Beyond FFT size and parameter dimension issues similar 
to item 1 and 2, particle filter design in Q needs to 
repeat the Kalman operation at multiple particle samples, 
yielding a multiplicative effect on complexity. 
To summarize, the major challenge in implementing a low- 
complexity EKF design lies in the factors of: (i) multiple 
TX users; (ii) multiple RX antennas; (iii) high parameter 
dimension; and (iv) large FFT size (e.g., 2048 size). 

With low-complexity and low-delay requirement in mind, 
this paper will present a time-sharing and interference miti- 
gation based Extended Kalman Filter (EKF) design for multi- 
user cognitive communications, which can estimate the CFO 
and channel parameters from multiple TX users to multiple RX 
antennas in a time-sharing manner. Here low delay property is 
achieved by using Kalman filter estimation at each RX sample 
in a real-time manner, and low complexity property is achieved 
by reusing a single user EKF design in a time-sharinjH and 
pipeling way. We first present a fundamental EKF design that 
estimates the CFO and channel parameters from one TX user 
to one RX antenna. Then such basic EKF design is reused in a 
time-shared way to estimate the parameters from multiple TX 
users to multiple RX antennas. Meanwhile, at each RX sample, 
an interference mitigation strategy is developed to estimate and 
remove the expected multi-user interference. In addition, we 
provide an adaptive noise variance tracking module to further 
enhance the estimation performance. Because of the usage 
of EKF structure, our design is essentially different from the 
particle filters in Q and the EM method in J5). Our design 
is also different from the Parallel-EKF design in [7| and the 
FFT-Block EKF design in O at the following perspectives: 
(i) our design runs at each time domain RX sample, not at an 
FFT-block basis; (ii) our design treats each user separately, not 
jointly; (iii) our design can be implemented in a time-sharing 
way, which is less considered in [6| and Q; (iv) system 
model in our design is different from the ones in (6), Q by 
integrating CFO parameter into channel response (see Eqn. 
<(5j in section II). Analysis and simulations results validate 
that our proposed design can closely approach the Cramer- 
Rao bound, and has lower computational complexity than the 
ones in J6), JT). Finally, although cognitive communication is 
a typical application scenario for our proposed design, it is 
also applicable in many other multi-user systems that satisfy 
the conditions presented in section II. 

II. System Model 

A. Problem Formulation 

We consider a total of Q TX users in the cognitive network. 
One of them is the primary TX user (i.e., base station), 
and the rest are all secondary TX users. Primary TX user's 
transmission is based on a time division MAC protocol, where 
time is divided into different time frames with equal duration. 

2 Time-sharing in this paper indicates that the same hardware module can 
be reused by different processes at separate time slots. 



Secondary users can maintain time synchronization with the 
primary TX user by learning and synchronizing with its time 
frames. Each secondary RX user is equipped with Na multiple 
antennas to decode the packets. Without loss of generality, 
we assume that every TX user has only one spatial stream, 
and there exists Q < Na- Also, every TX user has a distinct 
training symboQ s q (n) with 1 < q < Q and < n < Np — 1. 
Here Np is the FFT size of OFDM system. 

Each TX user has an independent carrier frequency offset 
(CFO) that is caused by both the oscillator offset and the 
doppler offset. Denote TX user g's CFO value as e q . For 
a given secondary RX user, the channel from TX user q to 
the mth RX antenna of this secondary user is denoted as 
hq t m{p q i ' m ), V 1 < I < L max . Here L max is the number of time 
domain paths in the channel response, and pf m is an integer 
value representing the relative delay of the Zth path. We assume 
that all pf m values have already been determined at an early 
stage (e.g., using PN sequences at coarse synchronization). 
The received signal at the mth RX antenna is derived as: 

Q ( 2ne n\ Lmax 
Vm{n) = J^exp ij q J 2J h q . m {pf m )s q [(n - pf' m ) N F ] 

q=l ^ F ' 1=1 

+z m (n),0<n<N F -l (1) 

where (n — p1' m )N F = {( n ~ Pi' m ) mod Nf} is circular 
shift, and z m (n) is the background noise at the rath RX 
antenna. 

The task in this work is to estimate CFO parameter e q and 
channel parameter h q ^ m (pf m ) for all users (1 < q < Q) and 
all antennas (1 < m < Na)- Obviously, the optimal estimation 
is the solution to this maximum likelihood (ML) problem: 

N A 

min ^ \Vm(n)- 

X^exp [3 27 Z qU ) K™(pf m ) s q i( n ~Pt m )N F ] | 2 (2) 

q=l ^ F ' 1 = 1 

where e q and h q , m (pf m ) represent the estimated values. There 
are a total of (L meix NA + 1)Q parameters in Eqn. (O, which 
constitutes a high-dimensional parameter estimation problem. 

B. State-Space Formulation 

ML solution can generally approach the optimal perfor- 
mance but it requires huge computations, which are highly 
undesirable in most systems. Instead, this paper proposes 
an EKF design for the estimation of the CFO and channel 
parameters, which can sequentially update the estimation 
results at each RX sample, resulting in low buffer size and 
low estimation delay. Initially, it is straightforward to directly 
apply an EKF design at Eqn. (f2]) by building all CFO and 
channel parameters into one state vector, whose dimension is 
as high as {L^^Na + This method will significantly 
increase the complexity of the derived Kalman filter. With 

3 This unique training symbol can be determined according to either the 
unique user ID in the network, or the access order in the current time frame. 



such complexity consideration in mind, we first propose a low- 
dimensional EKF design that can estimate the parameters from 
one TX user to one RX antenna, which has only (L max + 1) 
parameters. Then we reuse this fundamental EKF design 
in a time-shared manner to estimate the parameters from 
multiple TX users to multiple RX antennas. In this way, high- 
dimensional parameters are estimated by sequentially reusing 
a low-dimensional estimator, which reduces the complexity of 
the proposed EKF design. 

We first present an RX signal formulation from the perspec- 
tive of TX user q and the mth RX antenna as: 



( ,2-K£ q n 
Vq, m \n) = exp I J — 



1=1 



h qiTn (pf m )s q \{n-pf m ) 



(3) 



here y 9jm (n) is extracted from y m (n) with the aid of interfer- 
ence mitigation module, and z q:m (n) represents the residual 
noise at TX user q and the mth RX antenna, which includes 
both the residual co-channel interference from other TX users 
and the background noise at the mth RX antenna. Additionally, 
the initial value of y q , m (n) without any interference mitigation 
is set to Uq^min) — y m (n). Details about the interference 
mitigation module will be given in section III. 
Now we define the associated state vector as: 



Xg, m (n) = [£„H g , m (n)] , 
. 2imE n 



H q ,m(n) = exp j 



(4) 

[hq, m (pl' m ), .-, ^g,m(pI'™ M )j5) 



The state-space model for TX user q and the mth RX antenna 
can be derived from (O as: 



X 9>m (rc) = /{X g>m (n - 1)} = Dg >m (£g)Xg, m (n - 1), (6) 
1 lxLmax 
Oi m „xi exp (j2ne q /N F ) lL max xL, 
Vq,m{n) = G 9>m (n)Xg !m (n) + z q>m (n), 
G g ,m(n) = [0, s q [(n - p 9 1 ' m )N F ] , s q [(n - P2' m )jv F 



D 
(8) 



[( n_ Pw)^ 



(9) 



Finally, it is worth mentioning that the derived state-space 
formulation is a nonlinear model, since there is a nonlinear 
component exp (j2n£q/Np) in the matrix ~D q _ m (e q ). 

III. Interference Mitigation based EKF Design 

This section sequentially describes (i) the basic EKF design 
that aims at only one TX user and one RX antenna, (ii) 
the interference mitigation module that cancels co-channel 
interference at each RX sample, (iii) the proposed adaptive 
noise variance tracking module, and (iv) the overall paradigm 
of the proposed design. 

A. Fundamental EKF Design 

The key idea behind the EKF design is using Jacobian 
derivation to linearize the nonlinear matrix D ? m (e g ) at local 



estimates: 

F qtm (n - 1) 



df (X,, ro (n-1)) 



dX g , m (n - 1) 
1 



,(n-l|n— 1) 
OlxL„ 



a(n - l)H^ m (n - l|n - 1) exp (a(n - 1)) I,L m ., x xL n 



a{n ~ 1) = j2TTSq tm (n - l\n - l)/N F 

where X g , m (rt — l\n — 1) represents the estimated state vec- 
tor after processing the (n — l)th RX sample. Based on dTob . 
the prediction steps in our fundamental EKF design are: 



(10) 
(11) 



Xg, m (n|n - 1) = 

~D q , m (£q.rn{n - l\n - l))X^ >m (n - l\n - 1), 
P q , m (n\n - 1) = 

Fg im (n - l)P,, m (n - l\n- l)F* TO (n - 1). 
And the updating steps are as follows: 



Kq tm (n) 



q,77l 



(n\n-l)Gf im (n)x 



[Gq, m {n)'P q , m {n\n - l)G* m {n) + <J 2 q m {n)] 1 



P g , m (n|n) = 

[I - K q>m (n)G q>m (n)] P g , m (n 
X gm (n|n) = X 9im (n|n - 1)+ 



K g>ro (n) y q , m (n) 



1 (n)X gjm (n|n - 1) 



(12) 
(13) 

(14) 
(15) 

(16) 



Here <r q m (n) represents the variance of the observation noise 
z q , m {n) in Eqn. (0). Also, CFO estimate £ q , m (n\n) in the state 
vector X 9im (n|n) should only use its real part as e 9jm (n|n) = 
Real [e g , m (n|n)]. Finally, the EKF design presented above is 
only used for the parameter estimation of one TX user and 
one RX antenna. This basic EKF design is then iterated in a 
time-shared manner to estimate the parameters of all TX users 
(1 < q < Q) and all RX antennas (1 < m < N A ). 

B. Interference Mitigation and Refined CFO Estimation 

Before describing the interference mitigation module, we 
first look at the refinement of the CFO estimates. Although 
TX user q has only one CFO parameter, our proposed EKF 
design can result in Na different estimates that are derived 
from Na RX antennas, which are denoted as ^ q: m{n\n — 1), 
1 < m < Na- As a result, we can use these N A different 
estimates to get a refined result e q (n\n— 1), which is calculated 
as: 



£ q {n\n - 1) = 

l/Pg,mM« _ 1 )(1,1) 



N A 



Er=i VlVNn - 1) (M) 



e q ,m{n\n - 1), (17) 



where P 9jm (n|n — l)(i,i) denotes P q , m (n\n — l)'s element 
located at the 1 st row and 1 st column. 

Recall that y q . m (n) involved in (fjjl and (fl6] l is extracted 
from the original RX signal y m {n) with the help of an 
interference mitigation strategy. Having derived the refined 
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Fig. 1. Block diagram of our proposed design. 



CFO estimate e q {n\n — 1), now the interference mitigation 
process can be applied at y q77n (n) as follows: 



y q ,m{n) = Um{n)- 

. 2irE: u (n\n — 1) 



exp j 



G„ !m (n)X Uim (n|n - 1). 

(18) 



C. Adaptive Noise Variance Tracking 

Since y q . m (n) is extracted from y m (n) via interference 
mitigation module, the variance of z q , m (n), (i.e., a 2 m (n) used 
in Eqn. (fj4"l>). is varying during the convergence process of the 
interference mitigation module, which has to be adaptively 
tracked. Such variance tracking is based on the following 
observation: 



E 



Uq,m{n) - G g , m (rc)Xg im (7i|n - 1) 



G q , m (n)P q , m (n\n - l)G£ m (n) + o* m (n). (19) 
Using (19), noise variance cr 2 , m (n) can be tracked as: 



1 - 6»+i 



— T ■ {max[e,, m (n),0]} 



e q , m (n) = y q , m {n) - G g , m (n)X 9jm (n|ri - 1) 
Gq,m(n)P q , m (n\n - i)G" m (n), 



(20) 



(21) 



where b — 0.99 is the decay factor used to exponentially 
weight the history values. 

D. Block Diagram 

The complete functional diagram of our proposed design 
is shown in Fig. Q] In this paradigm, received samples at all 
RX antennas are first processed in the interference mitigation 
module. Then the resultant samples are sequentially processed 
in the basic EKF module and noise variance tracking module. 
For the ease of description, all components in Fig. Q] are de- 
picted in a parallel manner. However, in practice these design 
components can be implemented in a time-shared manner, and 
only one single EKF module is physically required. 



TABLE I 

Complexity Comparison (Number of Complex Multiplications) 



Design Name 


Number of Complex Multiplications 


Proposed Design 


~ ^max + lO^max + l^max + 2 


T7nll Ctnt Q T^L^TT 

rull-Mate L,Ivr 


/n J r")>i ( T AT 1 \3 1 
U [Q OmaxJVA + J.J } 


rrl-rlloCK EJvr |o| 


o{N A N F {QN A y{L m ^ + iy 

+ (QA r A) 3 (imax + l) 3 } 


Parallel EKF Q 


o{Q*N F {L max + iy + QN*(L max + iy 

+QA r F (L max + l) 3 } 


Particle Filter 


0{N P Q'-'Np (I/max + 1)'' + NpQN F (L maM + l)' z 
+NpQN F (L mB , + l) 3 } 
Np is the number of particle samples. 


EM method (5) 


0{N L QN F L max + NlQNpL^J 
Nl is the number of iterations. 



IV. Simulations and Discussions 

A. Parameters Setup 

The OFDM system built in the simulation has a bandwidth 
of 20MHz and a FFT size Np = 2048. Each TX user's 
CFO value is independently and randomly generated within 
the range of [-2, 20 Wireless channels are generated using 
the SUI-3 channel model [8|, which has L max = 3 non-zero 
paths at the time domain. SNR in this paper is defined as 
the ratio of the signal power to the noise power at one RX 
antenna, i.e., SNR = o\ja 2 z where a 2 R is the total received 
signal power at one RX antenna that is coming from all 
TX users, and a 2 z is the power of the background noise. In 
the simulation, CFO and channel parameters are estimated 
using one OFDM training symbol with Np = 2048 samples. 
Estimation results are validated via the mean square error 
(MSE) performance. Specifically, MSE for channel estimation 
is defined as a normalized version: 

2 



MSE(ft, )m ) 



Ef=T 



(Pi )- h q,m(Pi ) 



£f=f 



,(pD 



(22) 



Cramer-Rao bounds for the MSE results of CFO and channel 
estimation can be derived according to (5) as: 

3Q 



CRB CF o(SNR) 

CRBchan(SNR) = 

B. Simulation Results 



2tt 2 • N F ■ SNR ■ N A ' 

Q 



2N F SNR' 



(23) 
(24) 



We consider a cognitive network with one primary link 
and three secondary links (a total of 4 links). We investigate 
the CFO and channel estimations at one secondary RX user 
with Na = 4 RX antennas. Without loss of generality, we 
assume that this secondary RX user has the same received 
power from all TX users. We plot the MSE results of CFO 
and channel estimation in Fig. |2] and Fig. [3] respectively. The 

4 Since integer frequency offsets can generally be estimated during the 
coarse synchronization stage, CFO value at fine synchronization stage is 
usually between -0.5 and 0.5. But here we use range 2 to demonstrate our 
design's estimation performance. 



results show that our proposed design can closely approach 
the Cramer-Rao bounds. In addition, we repeat our simulation 
by disabling interference mitigation module, or noise variance 
tracking module. The corresponding results (Fig. [2] and Fig. 
[3]) indicate that without interference mitigation, the estimation 
performance can be dramatically degraded. And without noise 
variance tracking module, there could be an error floor at 
high SNR values because of the inaccurate tracking of noise 
variance information. Using the values in Fig. [2] and Fig. [3] 
it is feasible to further investigate the BER/PER performance. 
But such discussions heavily depend on the designed receiver 
structure, which is omitted here for page limitation. 

C. Delay, Buffer Size and Complexity Analysis 

This subsection evaluates the issues of complexity, delay 
and buffer size in the considered designs. We first look at 
the complexity issue. In particular, we count the number of 
complex multiplications involved in our proposed EKF design, 
which is listed in TableU And for comparison, in that table, we 
also list the complexity results of Full-State EKF, FFT-Block 
EKF 0, Parallel EKF Q, Particle filter (7), and EM method 
0. Here Full-State EKF represents the EKF that builds all 
(L max NA + 1)Q states in (O into one state vector, yielding 
high state dimension. We can see that our proposed design 
enjoys the lowest computation complexity, which is only at 
the order of £f nax . But Full-State EKF's complexity is around 
Q^N\ higher than our design. Moreover, FFT-Block EKF, 
Parallel EKF, Particle Filter, and EM method's complexities 
all rely on FFT size Np, which is significantly large in our 
case (N F = 2048). 

Now we further look at the delay and buffer size in the 
proposed design. Since our EKF scheme updates the Kalman 
estimate at each RX sample (not at each FFT block) in an 
online and real-time manner, it has low estimation delay and 
requires low buffer size. However, Particle Filter [7|, Parallel 
EKF [7 1, and EM approach [5| all operate at an FFT-block 
basis with buffer size Np = 2048 samples, resulting in both a 
large delay and a large buffer size. Even worse, particle filter 
and EM method both need to process the FFT-block multiple 
times (e.g., particle samples in particle filter, and iterations in 
EM method), leading to additional estimation delay. 

V. Conclusion 

This paper has presented a low-delay and low-complexity 
EKF design that can estimate the CFO and channel parameters 
in multi-user cognitive communications. We first present a 
fundamental EKF design that works at one TX user and 
one RX antenna. Then this basic EKF design is reused in 
a time-shared way to estimate the parameters for multiple 
TX users at multiple RX antennas. Besides, an interference 
mitigation strategy is proposed to estimate and cancel the 
multi-user interference at each RX sample. Moreover, adaptive 
noise variance tracking module is further employed to further 
enhance the estimation performance. Compared with existing 

5 Particle filters in |7 | and EM designs in (5) have even higher complexity 
because of either the number of particle samples, or the number of iterations. 
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Fig. 2. CFO estimation's MSE results under different SNR values. 
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Fig. 3. Channel estimation's MSE results under different SNR values. 



related designs (FFT-Block EKF 0, Parallel EKF Q, Particle 
filter [7 1, and EM method O), our proposed design enjoys 
low computation complexity (because of pipelining and time- 
sharing design), low delay and low buffer size (due to its online 
and run-time estimation). Besides, its estimation performance 
can closely approach the Cramer-Rao bound. 
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Abstract — Parameter estimation in cognitive communications 
can be formulated as a multi-user estimation problem, which is 
solvable under maximum likelihood solution but involves high 
computational complexity. This paper presents a time-sharing 
and interference mitigation based EKF (Extended Kalman Filter) 
design for joint CFO (carrier frequency offset) and channel 
estimation at multiple cognitive users. The key objective is to 
realize low implementation complexity by decomposing high- 
dimensional parameters into multiple separate low-dimensional 
estimation problems, which can be solved in a time-shared 
manner via pipelining operation. We first present a basic EKF 
design that estimates the parameters from one TX user to one 
RX antenna. Then such basic design is time-shared and reused 
to estimate parameters from multiple TX users to multiple RX 
antennas. Meanwhile, we use interference mitigation module 
to cancel the co-channel interference at each RX sample. In 
addition, we further propose adaptive noise variance tracking 
module to improve the estimation performance. The proposed 
design enjoys low delay and low buffer size (because of its online 
real-time processing), as well as low implementation complexity 
(because of time-sharing and pipeling design). Its estimation 
performance is verified to be close to Cramer-Rao bound. 

I. Introduction 

Cognitive communication system is widely accepted as 
a perspective way in increasing the spectrum efficiency of 
wireless networks, where primary links and secondary links 
can usually co-exist in the network, resulting in an interference 
limited environment. Parameter estimation in cognitive com- 
munications is a challenging problem because of (i) the exis- 
tence of co-channel interference, and (ii) the high-dimensional 
parameters from multiple TX users to multiple RX anten- 
nas. In particular, note that different TX users often have 
independent carrier frequency offset (CFO) values (including 
both oscillator offsets and doppler offsets), which usually 
introduce serious nonlinear components within the observed 
signal, complicating the estimation problem. Meanwhile, chan- 
nel responses from multiple TX users to multiple RX antennas 
can result in a set of high-dimensional parameters, which are 
also difficult to estimate. Finally, due to the existence of multi- 
user interference, CFO and channel parameters usually have 
to be treated together and be estimated in a joint way so as 
to approach the optimal performance, which further increases 
the estimation complexity. 

Without loss of generality, this paper assumes Orthogonal 
Frequency Division Multiplexing (OFDM) system, which is 



an overwhelming choice for modern wireless systems. The 
classical CFO and channel estimation method in a single-user 
OFDM system is based on two repeated training symbols 
[?]. It has low implementation complexity and near-optimal 
performance, but only applies to a single-user scenaricQ. In 
multi-user OFDMA systems with unique subcarrier set per 
TX user, CFO and channel parameters can be recovered by 
exploiting distinct subcarrier structures among TX users ([?], 
[?]). But this method requires separate subcarrier allocation for 
different users. Consequently, in a general multi-user cognitive 
system without specific subcarrier allocation per user and with 
overlapped training symbols, ML and EM related methods 
seem to be the only applicable choice, where all TX users' 
parameters have to be formulated into a maximum likelihood 
(ML) estimation problem [?], which is solvable under Ex- 
pectation Maximization (EM) method [?] in an iterative way. 
However, such ML and EM approaches often require high 
computational complexity and high processing delay, because 
in these methods, the entire OFDM block is stored offline and 
is iteratively processed multiple times. 

Based on above considerations, this paper will focus on 
using Kalman filter structure to estimate the CFO and channel 
parameters in multi-user cognitive communications. Our major 
objective is to achieve low-complexity and low-delay estima- 
tion performance in cognitive systems. In general, Kalman 
filter is a good candidate for low delay and low complexity pa- 
rameter estimation primarily due to its real-timing processing 
property. It has been conventionally used for CFO and channel 
estimation in multi-user OFDM systems, e.g., the FFT-Block 
EKF design in [?], the parallel EKF design in [?], and the 
particle filter design in [?]. However, these existing designs 
inherently suffer from multiple issues related to complexity, 
delay and buffer size as follows: 

1) Block EKF design in [?] operates on an FFT-block 
basis, which grows increasingly complex as FFT size 
becomes large (e.g., 2048 FFT size). Also, parameters 
in this method are jointly estimated with high parameter 
dimension. 

2) Parallel EKF design in [?] also operates on an FFT-block 
basis, marking it complex under large FFT size. Parame- 

1 It is also applicable for multi-user scenario with non-overlapping training 
symbols, but this is not the case considered in this paper. 



ters in this design are jointly estimated by calculating the 
covariance information between different users, leading 
to a high matrix dimension. 
3) Beyond FFT size and parameter dimension issues similar 
to item 1 and 2, particle filter design in [?] needs to 
repeat the Kalman operation at multiple particle samples, 
yielding a multiplicative effect on complexity. 
To summarize, the major challenge in implementing a low- 
complexity EKF design lies in the factors of: (i) multiple 
TX users; (ii) multiple RX antennas; (iii) high parameter 
dimension; and (iv) large FFT size (e.g., 2048 size). 

With low-complexity and low-delay requirement in mind, 
this paper will present a time-sharing and interference miti- 
gation based Extended Kalman Filter (EKF) design for multi- 
user cognitive communications, which can estimate the CFO 
and channel parameters from multiple TX users to multiple RX 
antennas in a time-sharing^manner. Here low delay property is 
achieved by using Kalman filter estimation at each RX sample 
in a real-time manner, and low complexity property is achieved 
by reusing a single user EKF design in a time-sharing and 
pipeling way. We first present a fundamental EKF design that 
estimates the CFO and channel parameters from one TX user 
to one RX antenna. Then such basic EKF design is reused in a 
time-shared way to estimate the parameters from multiple TX 
users to multiple RX antennas. Meanwhile, at each RX sample, 
an interference mitigation strategy is developed to estimate and 
remove the expected multi-user interference. In addition, we 
provide an adaptive noise variance tracking module to further 
enhance the estimation performance. Because of the usage 
of EKF structure, our design is essentially different from the 
particle filters in [?] and the EM method in [?]. Our design 
is also different from the Parallel-EKF design in [?] and the 
FFT-Block EKF design in [?] at the following perspectives: 
(i) our design runs at each time domain RX sample, not at an 
FFT-block basis; (ii) our design treats each user separately, not 
jointly; (iii) our design can be implemented in a time-sharing 
way, which is less considered in [?] and [?]; (iv) system 
model in our design is different from the ones in [?], [?] by 
integrating CFO parameter into channel response (see Eqn. 
<(5j in section II). Analysis and simulations results validate 
that our proposed design can closely approach the Cramer- 
Rao bound, and has lower computational complexity than the 
ones in [?], [?]. Finally, although cognitive communication is 
a typical application scenario for our proposed design, it is 
also applicable in many other multi-user systems that satisfy 
the conditions presented in section II. 

II. System Model 

A. Problem Formulation 

We consider a total of Q TX users in the cognitive network. 
One of them is the primary TX user (i.e., base station), 
and the rest are all secondary TX users. Primary TX user's 
transmission is based on a time division MAC protocol, where 

2 Time-sharing in this paper indicates that the same hardware module can 
be reused for different processes at separate time slots. 



time is divided into different time frames with equal duration. 
Secondary users can maintain time synchronization with the 
primary TX user by learning and synchronizing with its time 
frames. Each secondary RX user is equipped with Na multiple 
antennas to decode the packets. Without loss of generality, 
we assume that every TX user has only one spatial stream, 
and there exists Q < Na- Also, every TX user has a distinct 
training symboQ s q (n) with 1 < q < Q and < n < Np — 1. 
Here Np is the FFT size of OFDM system. 

Each TX user has an independent carrier frequency offset 
(CFO) that is caused by both the oscillator offset and the 
doppler offset. Denote TX user q's CFO value as e q . For 
a given secondary RX user, the channel from TX user q to 
the mth RX antenna of this secondary user is denoted as 
h q ,m(pf m ), V 1 < I < L max . Here i max is the number of time 
domain paths in the channel response, and pj' m is an integer 
value representing the relative delay of the Zth path. We assume 
that all pf' m values have already been determined at an early 
stage (e.g., using PN sequences at coarse synchronization). 
The received signal at the mth RX antenna is derived as: 

Q / 2ire n\ Lmax 

Vm{n) = ^ eX P U^r - ) X! h Q,m(Pl' m )Sq [{n ~ pf^Np] 
9=1 V F ' 1=1 

+z m {n),0 < n < N F - 1 (1) 

where (n — P 9 ' m )N F = {{n — Pi' m ) mod Np} is circular 
shift, and z m (n) is the background noise at the mth RX 
antenna. 

The task in this work is to estimate CFO parameter e q and 
channel parameter h q , m (p1' m ) for all users (1 < q < Q) and 
all antennas (1 < m < Na)- Obviously, the optimal estimation 
is the solution to this maximum likelihood (ML) problem: 

N A 

mill ^ \Vm{n)- 

m— 1 

X^exp ( 3 27 ^ qU ) K™{pT m )s q [{n-pf m ) NF ] | 2 (2) 

9=1 ^ F ' 1 = 1 

where e q and h qt7n (p'j l ' m ) represent the estimated values. There 
are a total of (L meix NA + 1)Q parameters in Eqn. (0, which 
constitutes a high-dimensional parameter estimation problem. 

B. State-Space Formulation 

ML solution can generally approach the optimal perfor- 
mance but it requires huge computations, which are highly 
undesirable in most systems. Instead, this paper proposes 
an EKF design for the estimation of the CFO and channel 
parameters, which can sequentially update the estimation 
results at each RX sample, resulting in low buffer size and 
low estimation delay. Initially, it is straightforward to directly 
apply an EKF design at Eqn. (fJI by building all CFO and 
channel parameters into one state vector, whose dimension is 
as high as (L max NA + This method will significantly 

3 This unique training symbol can be determined according to either the 
unique user ID in the network, or the access order in the current time frame. 



increase the complexity of the derived Kalman filter. With 
such complexity consideration in mind, we first propose a low- 
dimensional EKF design that can estimate the parameters from 
one TX user to one RX antenna, which has only (L max + 1) 
parameters. Then we reuse this fundamental EKF design 
in a time-shared manner to estimate the parameters from 
multiple TX users to multiple RX antennas. In this way, high- 
dimensional parameters are estimated by sequentially reusing 
a low-dimensional estimator, which reduces the complexity of 
the proposed EKF design. 

We first present an RX signal formulation from the perspec- 
tive of TX user q and the mth RX antenna as: 
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here y 9jm (n) is extracted from y m (n) with the aid of interfer- 
ence mitigation module, and z q:m (n) represents the residual 
noise at TX user q and the mth RX antenna, which includes 
both the residual co-channel interference from other TX users 
and the background noise at the mth RX antenna. Additionally, 
the initial value of y q ,m(n) without any interference mitigation 
is set to y q ,m{n) — y m (n). Details about the interference 
mitigation module will be given in section III. 
Now we define the associated state vector as: 
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The state-space model for TX user q and the mth RX antenna 
can be derived from (O as: 
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Finally, it is worth mentioning that the derived state-space 
formulation is a nonlinear model, since there is a nonlinear 
component exp (j2Tre q /N F ) in the matrix T> q ^ m {e q ). 

III. Interference Mitigation based EKF Design 

This section sequentially describes (i) the basic EKF design 
that aims at only one TX user and one RX antenna, (ii) 
the interference mitigation module that cancels co-channel 
interference at each RX sample, (iii) the proposed adaptive 
noise variance tracking module, and (iv) the overall paradigm 
of the proposed design. 

A. Fundamental EKF Design 

The key idea behind the EKF design is using Jacobian 
derivation to linearize the nonlinear matrix D ? m (e g ) at local 
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where X g , m (rt — l|n — 1) represents the estimated state vec- 
tor after processing the (n — l)th RX sample. Based on dlOl ). 
the prediction steps in our fundamental EKF design are: 
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P q , m (n\n - 1) = 

Fg im (n - l)P,, m (n - l\n- l)F* TO (n - 1). 
And the updating steps are as follows: 
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Here <r q m (n) represents the variance of the observation noise 
z q , m {n) in Eqn. (0). Also, CFO estimate e q , m (n\n) in the state 
vector X 9im (n|n) should only use its real part as e 9jm (n|n) = 
Real [e g , m (n|n)]. Finally, the EKF design presented above is 
only used for the parameter estimation of one TX user and 
one RX antenna. This basic EKF design is then iterated in a 
time-shared manner to estimate the parameters of all TX users 
(1 < q < Q) and all RX antennas (1 < m < N A ). 

B. Interference Mitigation and Refined CFO Estimation 

Before describing the interference mitigation module, we 
first look at the refinement of the CFO estimates. Although 
TX user q has only one CFO parameter, our proposed EKF 
design can result in Na different estimates that are derived 
from Na RX antennas, which are denoted as ^ q: m{n\n — 1), 
1 < m < Na- As a result, we can use these N A different 
estimates to get a refined result e q (n\n— 1), which is calculated 
as: 
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e q , m {n\n - 1), (17) 



where P 9jm (n|n — l)(i,i) denotes P q , m (n\n — l)'s element 
located at the 1 st row and 1 st column. 

Recall that y q . m (n) involved in (0 and (fTSI l is extracted 
from the original RX signal y m {n) with the help of an 
interference mitigation strategy. Having derived the refined 
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Fig. 1. Block diagram of our proposed design. 



CFO estimate e q {n\n — 1), now the interference mitigation 
process can be applied at y q77n (n) as follows: 



y q ,m{n) = Um{n)- 

. 2irE: u (n\n — 1) 



exp j 
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C. Adaptive Noise Variance Tracking 

Since y q . m (n) is extracted from y m (n) via interference 
mitigation module, the variance of z q , m (n), (i.e., a 2 m (n) used 
in Eqn. (TT4"1>). is varying during the convergence process of the 
interference mitigation module, which has to be adaptively 
tracked. Such variance tracking is based on the following 
observation: 
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Using (19), noise variance cr 2 , m (n) can be tracked as: 
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where b — 0.99 is the decay factor used to exponentially 
weight the history values. 

D. Block Diagram 

The complete functional diagram of our proposed design 
is shown in Fig. Q] In this paradigm, received samples at all 
RX antennas are first processed in the interference mitigation 
module. Then the resultant samples are sequentially processed 
in the basic EKF module and noise variance tracking module. 
For the ease of description, all components in Fig. Q] are de- 
picted in a parallel manner. However, in practice these design 
components can be implemented in a time-shared manner, and 
only one single EKF module is physically required. 



TABLE I 

Complexity Comparison (Number of Complex Multiplications) 



Design Name 


Number of Complex Multiplications 
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Parallel EKF [?] 
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Particle Filter [?] 
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Np is the number of particle samples. 


EM method [?] 
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Nl is the number of iterations. 



IV. Simulations and Discussions 

A. Parameters Setup 

The OFDM system built in the simulation has a bandwidth 
of 20MHz and a FFT size Np = 2048. Each TX user's 
CFO value is independently and randomly generated within 
the range of [-2, 20 Wireless channels are generated using 
the SUI-3 channel model [?], which has L ma x = 3 non-zero 
paths at the time domain. SNR in this paper is defined as 
the ratio of the signal power to the noise power at one RX 
antenna, i.e., SNR = cr\/a 2 z where a 2 R is the total received 
signal power at one RX antenna that is coming from all 
TX users, and a 2 z is the power of the background noise. In 
the simulation, CFO and channel parameters are estimated 
using one OFDM training symbol with Np = 2048 samples. 
Estimation results are validated via the mean square error 
(MSE) performance. Specifically, MSE for channel estimation 
is defined as a normalized version: 
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Cramer-Rao bounds for the MSE results of CFO and channel 
estimation can be derived according to [?] as: 
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B. Simulation Results 



2tt 2 • N F ■ SNR ■ N A ' 

Q 



2N F SNR' 



(23) 
(24) 



We consider a cognitive network with one primary link 
and three secondary links (a total of 4 links). We investigate 
the CFO and channel estimations at one secondary RX user 
with Na = 4 RX antennas. Without loss of generality, we 
assume that this secondary RX user has the same received 
power from all TX users. We plot the MSE results of CFO 
and channel estimation in Fig. |2] and Fig. [3] respectively. The 

4 CFO value at fine synchronization stage is generally between -0.5 and 
0.5, because any integer frequency offset can be estimated during the coarse 
synchronization stage. But here we use range 2 to demonstrate our design's 
estimation performance. 



results show that our proposed design can closely approach 
the Cramer-Rao bounds. In addition, we repeat our simulation 
by disabling interference mitigation module, or noise variance 
tracking module. The corresponding results (Fig. [2] and Fig. 
[3]) indicate that without interference mitigation, the estimation 
performance can be dramatically degraded. And without noise 
variance tracking module, there could be an error floor at 
high SNR values because of the inaccurate tracking of noise 
variance information. Using the values in Fig. [2] and Fig. [3] 
we can further investigate the BER/PER performance in the 
system. But such discussions heavily depend on the designed 
receiver structure, which is omitted here for page limitation. 

C. Delay, Buffer Size and Complexity Analysis 

This subsection evaluates the issues of complexity, delay 
and buffer size in the considered designs. We first look at 
the complexity issue. In particular, we count the number of 
complex multiplications involved in our proposed EKF design, 
which is listed in TableU And for comparison, in that table, we 
also list the complexity results of Full-State EKF, FFT-Block 
EKF [?], Parallel EKF [?], Particle filter [?], and EM method 
[?]. Here Full-State EKF represents the EKF that builds all 
{L max NA + 1)Q states in (O into one state vector, yielding 
high state dimension. We can see that our proposed design 
enjoys the lowest computation complexity, which is only at 
the order of £f nax . But Full-State EKF's complexity is around 
Q Z N\ higher than our design. Moreover, FFT-Block EKF, 
Parallel EKF, Particle Filter, and EM method's complexities 
all rely on FFT size Np, which is significantly large in our 
case (Np = 2048). 

Now we further look at the delay and buffer size in the 
proposed design. Since our EKF scheme updates the Kalman 
estimate at each RX sample (not at each FFT block) in an 
online and real-time manner, it has low estimation delay and 
requires low buffer size. However, Particle Filter [?], Parallel 
EKF [?], and EM approach [?] all operate at an FFT-block 
basis with buffer size Np — 2048 samples, resulting in both a 
large delay and a large buffer size. Even worse, particle filter 
and EM method both need to process the FFT-block multiple 
times (e.g., particle samples in particle filter, and iterations in 
EM method), leading to additional estimation delay. 

V. Conclusion 

This paper has presented a low-delay and low-complexity 
EKF design that can estimate the CFO and channel parameters 
in multi-user cognitive communications. We first present a 
fundamental EKF design that works at one TX user and one 
RX antenna. Then this basic EKF design is reused in a time- 
shared way to estimate the parameters for multiple TX users 
at multiple RX antennas. Besides, an interference mitigation 
strategy is proposed to estimate and cancel the multi-user 
interference at each RX sample. Moreover, adaptive noise 
variance tracking module is also employed to further enhance 
the estimation performance. Compared with existing related 
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Fig. 2. CFO estimation's MSE results under different SNR values. 
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Fig. 3. Channel estimation's MSE results under different SNR values. 



designs (FFT-Block EKF [?], Parallel EKF [?], Particle filter 
[?], and EM method [?]), our design enjoys low computation 
complexity (because of pipelining and time-sharing design), 
as well as low delay and low buffer size (due to its online and 
run-time estimation). Meanwhile, its estimation performance 
can closely approach the Cramer-Rao bound. 



5 Particle filters in [?] and EM designs in [?] have even higher complexity 
because of either the number of particle samples, or the number of iterations. 



